THE GEOMETRY OF SAMPLING ON UNIONS OF LATTICES 
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Abstract. In this short note we show two results concerning sampHng translation invariant 
subspaces of I/^(R'^) on unions of lattices. The first result shows that the sampling transform 
on a union of lattices is a constant times an isometry if and only if the sampling transform on 
each individual lattice is so. The second result demonstrates that the sampling transforms 
of two unions of lattices on two bands have orthogonal ranges if and only if correspondingly 
the sampling transforms of each pair of lattices have orthogonal ranges. We then consider 
sampling on shifted lattices. 



1. Introduction 

In this note we consider sampling translation invariant subspaces on unions of lattices in 
L^(]R'^). The translation invariant subspaces have the form 

Ve := {/ G L2(M'^) : suppU) C E} 

where £^ is a band, i.e. a measurable subset of finite measure. We normalize the Fourier 
transform so that for / G L^{W^) H L'^{W^), 



ho = / f{x)e-'-'^<dx. 



Define, for any y G M.'^, the translation operator Tyf{x) = f{x — y). Note that TyVE = Ve for 
all y and all E. Since Ve is a reproducing kernel Hilbert space, we have that, via the inverse 
Fourier transform, 

(1) /(y) = (/,T,^E), 

where ipE = Xe, the indicator function of the set E. Let {Ai, . . . , An} he dxd invertible (real) 
matrices. Define the sampling transform of the samples A := {AjZ+Pj : j = 1, . . . ,n; z E Z"^} 

by 

Qa-.Ve^ ©^=i/'(Z'^) : / ^ {fiA.z + A), . . . , f{AnZ + 

provided 6^ is bounded. A straight forward computation (see P) shows that 6a is bounded 
if and only if for each j = 1, . . . , n. 

We say the samples {AjZ + j3j : j = 1, . . . ,n] z G Z*^} is a set of sampling for the band 
E if 0A : Ve — > ©"^^^/^(Z'^) as defined above is bounded above and below. The samples 
{AjZ + : J = l,...,n; z e Z^} on the band E are tight if He^/f = K\\f\\^ for all 
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f G Ve- The samples are an exact set of sampling for the band E if the range of is all of 

If we oversample the space Ve, then the range of the sampling transform 0^^ will have 
a non-trivial orthogonal complement in t^il/) © • ■ ■ © t^ilJ^). In several applications, one 
wishes to know the range of the sampling transform. For example, this information is useful 
for denoising applications (see Additionally, for multiple access communications, one also 
wishes to know the range of the sampling transform (see [Ij). We address this in theorems 
2,3 and 4, by characterizing when two sampling transforms have orthogonal ranges. 

Definition 1. Suppose E and F are hands in W^. Let Ai, . . . An and Bi, . . . , Bn be invertible 
dx d matrices, and . . . , /3„, 71, ... , 7„} C M"'. Let A := {Ajz + Pj : z E Z'^, j = 1, . . . ,n} 
and B := {Bjz + 'jj : z E Z'*, j = 1, . . . ,n} be such that the sampling transforms Qa o-nd 
B o,re bounded on E and F , respectively. We say the samples A and B are orthogonal on 
E and F if QaVe is orthogonal to 0bVf ffi"=i/^(Z'^); equivalently 0^0^ = 0. 

We will use the following notation. If C is an invertible matrix, denote C = C*^^, where 
C* is the (conjugate) transpose. If f,g E L^(R°'), define the C bracket product to be 

[f,9]c{0 = E f{^ + Cm)g{^ + Cm). 

Let n'^ denote the unit cube [0, l]'^ C M'^. Finally, define 

V = {f e L2(M^) : / G L°°(M'^) and compactly supported}; 
clearly D is a dense subspace of ^^(R'^). 

2. Unshifted Lattices 

We set out to prove the following two theorems, the main results of the paper, which relate 
the behavior of sampling on a union of (unshifted) lattices to the behavior of sampling on 
individual lattices. 

Theorem 1. Let E be a band in M'^, Ai, . . . ,An be d x d matrices, and suppose that the 
sampling transform 0a is bounded on Ve- The samples {Ajz : j = 1, . . . , n; Z*^} are tight on 
E with constant K if and only if for each j the samples {AjZ : z G Z"^} are tight on E, with 
constant Kj. In this case, we have Kj = \ det Aj\~^ and K = X]j=i ^j- 

Theorem 2. Let E and F be bands in Mf^, Ai, . . . ,An and Bi, . . . ,Bn be d x d invertible 
matrices, and suppose that Qa o-nd Qb are bounded on Ve and Vp, respectively. The samples 
{AjZ : j = 1, . . . ,n; z E Z'^} and {BjZ : j = 1, . . . ,n] z E Z*^} are orthogonal on the bands E 
and F if and only if for j = 1, . . . , r 

(2) Yl + P)) • E ^^(^K^ + ?)) = a.e. C 

Equivalently, the samples above are orthogonal if and only if the samples {Ajz : z G Z'^} and 
{BjZ : z G Z'^} are orthogonal on the bands E and F for each j = 1, . . . ,r. 

As in jlj, let V he a. countable index set, let Cphe a. d x d invertible matrix for each p E V, 
and define the following: 

A = UpevC'X 
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and for a G A, 

V^ = {peV: C;a G Z^}. 
Note that if a = C'p^z for some z E \ {0}, then pq G Va', if a = 0, then Va = V. Let 
{gp:peV}c L'^iW^). The collection {Tc^kQp : p G P, A; G Z'^} is said to be Bessel if there 
exists a constant M < oo such that for all / G //^(M*^), 

pe-p k&<^ 

The collection {Tcpkdp : p E V, k E Z"^} satisfies the local boundedness condition if for every 
^(/)^=EE / MX^ + C',k)\'\detC,\-'\g,iO\'d^<^- 



Lemma 1. Suppose the sets G := {TcpkQp : p E V, k E Z*^} and H := {Tc^khp : p E V, k E 
Z'^} are both Bessel and both satisfy the local boundedness condition L{f) < oo for all f eT), 
as defined above. Define the operator 

nH,G : L\R') ^ L'iR") : / ^ ^ (/, Tcpk9p)Tc,khp. 

per feezd 

The operator Qh,g commutes with the translation operators Ty for all y eW^ if and only if 
for alia E A \ {0}, 

^ I det Cp\'^gp{^)hp{^ + a) = a.e. 
peVc 

In this case, Qh.g is a Fourier multiplier with symbol 

.(0:=5^|detCp|-i^/z,(0, 

per 

I.e. n^f{0 = s{OhO- 

Proof. For f E V, define the continuous function 

Wf{x) = {VtH^cT^f.T^f). 

If VLh,g commutes with all for x E M'^, then clearly Wf{x) is constant for all f E V. 
Conversely, if Wf{x) is constant for all f E V, then {T_xflH,GTxf, f) = (^_h",g/)/)) whence 
by the polarization identity, T_x^h,gTx = ^h,g, and thus VLh,gTx = Tx^h,g- 

By ^ Proposition 2.4], Wf{x) coincides pointwise with the almost periodic function 

where 

[ + V I det Cp\-'^hp{^ + a)di. 



Wf[a) 



By 131 Lemma 2.5] and the proof of Theorem 2.1 in |4^, Wf{x) is constant for all / G P if and 
only if for all a G A \ {0}, 

J2 I det Cp\-^'^hp{^ + a) = Oa.e. ^ 
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It is well known that if ^Ih,g commutes with Ty for all y G M*^, then it is a Fourier multiplier. 
Evaluating Wf{x) ai x = yields 

J^" pep 

Therefore, the symbol of ^Ih,g is s{^) as above. □ 
Lemma 2. Suppose {gp}pev C and let G be as in lemma{J[ Define Qq 

Qg ■■ L\R'') ^ xV):f^ {{f,Tc^kgp))(k,p)- 

We have HOcfi'lP = -^IkP fof' o-^l 9 ^ Ve 'if only if for all a G A, 

(3) Yl I detCpr^^^p(e + a) = K5^,oXe{0 a.e. 

where 5^,0 ^-5 the Kronecker delta. 

Proof. Note that, as defined in lemmaHl ^g,g = ©g^G) hence ||0g5'|P = {^G,Gg,g)- 
Suppose equation El holds. Then, by lemma HJ fl^G is a Fourier multipher, with symbol 

3(0 = 5^ I detCp|-i|^p(OP = 
per 

Thus, 

liec^^f = {nG,G9,9) = {s{Om,m) = {Kg,g) = K\\gf. 

Conversely, suppose IIOGf?!!^ = -^llfi'll^ for all g G Ve- Then, by the above computation, 
it follows that {QG,Gg,g) = {Kg,g), whence by the polarization identity, flG,G is a Fourier 
multiplier with symbol s{^) = K a.e. C,- It now follows that equation El holds. □ 

We wish to apply lemmas [H and |21 to theorems [H and |21 Let A := {Ajz + Pj : j = 
1, . . . , n; 2 G Z'^}, where Aj is an invertible matrix for each j, and let be a band such that 
Qa is bounded on Ve. By equation [H the sampling transform Qa above can be written in 
terms of the collection 

(4) {TA^,Tp^iJE:j = l,...,n; zeZ"}, 
as follows: 

ifiA.z + A), . . . , f{A^z + P^)) = (if, Ta.^TpM, ...,(/, Ta^.Tp^^e)). 

Hence, questions of tightness and orthogonality of sampling can be stated in terms of the 
collection (jD). In order to apply the lemmas, we need to verify that (jH) satisfies the local 
boundedness condition. 

Lemma 3. Let A, E be as above. If the sampling transform Qa is bounded on Ve, then for 
all f we have 

n „ 

^(/) = EE / .\h^ + A',m)\'\detA,nxE{0\'d^<^- 



THE GEOMETRY OF SAMPLING ON UNIONS OF LATTICES 

Proof. Let Q = ||/||oo- Since Qa is bounded on Ve, then the cardinahty of the set 



7d 

is essentially bounded in ^. Let this bound be R. Clearly, we have 



{m G Z*^ : ^ + A'-m G E} 



n „ 

L{f) < I detAj\-^ / QRd^< oo. 
i=i 



□ 

The following lemma is a variation of a lemma in it will be needed for theorem |2l 

Lemma 4. Let G be as in lemma QJ- let Dp be a d x d invertible matrix for each p & V, 
and suppose H := {Tp^^hp : p E V, k E U^} is Bessel and satisfies the local boundedness 
condition. Define ^h,g analogous to lemma^- 

For all f.geV, 

Proof. We have via the Fourier transform 

{^H,Gf, 9) = YY1 ^Cpk9p) {TDpkhp, 9) 

Each of the above integrals can be written as follows, by variable substitutions. 
f ^"(^)g2.iC,K-;(^rf^=|detC,|-i / f{C'pOe'''"'-^'WiX)d^ 

= J2 \ detCp\-' [ f{C'^{i + m))e'-^'^-%{C'p{i + m))di 

= |detC,ri / {l9v\c',{C'pi)e^'^'''<di. 
The interchanging of the sum and integral is valid since f,gED. Therefore, we have 
(5) {^H,Gf,9) = J2 \AJc D\ ^ / ^f'9ph^^p0e'''"'-^d^- [ [hp, g]n'^iD'^co)e~'-^'-du 
where the sum over k is really the inner product of the Fourier coefficients of the two functions 



[f.9p]c'^{C'pi) and [hp,g\D'^{D'^uj). 



-'pK p^J L P'yjJ^pV^p 

Thus, the sum (jSj) is 



pev 
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□ 

Proof of Theorem 1. Let V = {1, • • • ,n}, Cj = Aj, and gj = ipE for j = 1, . . . ,n. Thus, 
A = u]^-^A'jZ'^ and for a e A, Va = {j e {1, . . . ,n} : A*a G Z'^} as above. Note that if 
a = A'j^^z for some z E Z'^ \ {0}, then jo ^ "Pa- By equation [T] and lemmas |21 and El we have 
that 

ii0^/ir = ^ii/ir 

if and only if 

^ \detAj\-\EiOXEi^ + a) = K6^,oXEiO 

for all a G A. 

It follows that for z ^Z^\ {0} and for each j G {1, . . . , n}, 

+ = Oa.e. i. 

Therefore, for each j, 

Y,XE{^,{i^ra))<\a.e. ^ 
Since the sum is bounded above by 1, for all / G Ve H P, 

iie^^./f = idetA,r^ii/ir, 

(see and [T]). □ 

Proof of Theorem ^. In we show that equation |21 is sufficient for QaVe -L QbVf; thus we 
wish to show that it is necessary. Let Eq G E and Fq G F he any measurable bounded sets. 
Define go G Ve and /o G Vf by go = xeq and /o = X-Fo- Suppose QaVe -L ©b^f; by lemma El 
we have 

n „ 

^|detA,i?,|-' / [xe,,Xe\a'^{^,0[xf,Xf,]b'^{B'^0=^- 
i=i 

Therefore, for each j = 1, . . . , r, 

/ [xeo,Xe]a'^{A'^0[Xf,Xfo]b'^{B',0 =0- 
However, expanding the above integrand yields 

I E Xi.oK(e + P)) Yl XFAB-{i + q)) = 0. 
The result now follows since Eo and Fq were arbitrary measurable bounded subsets. □ 



3. Shifted Lattices 

We begin with the following remark. Both theorems 1 and 2 are no longer valid if the 
lattices are allowed to be shifted. As a trivial example for theorem 1, consider the canonical 
bandlimited function space Vi_i/2,i/2] and the lattices 2Z and 2Z + 1; individually neither form 
sets of sampling, but together they form an exact set of sampling. An example corresponding 
to theorem 2 is as follows. 
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Example 1. Let E = [0, 1], F — [—1,0] and define the following sampling transforms on Ve 
and Vp respectively: 

Qa-.Ve^ © l\Z) : / ^ (/(n). f{n + 1/2)) 
Qb-Vf^ l\Z) © . g ^ (^(^ + 1/2), g{n)). 

Define 0yi^, and 0^^, for j = 1,2, to be the first and second coordinates of 0^ and 0^ 
respectively. It is clear that QajVe = QbjVf = /^(Z). However, we claim that the ranges of 
0^ and Qb are orthogonal. 

Define a unitary operator U : Ve ^ Vp : f{C) i— > f{C + 1)- By the polarization identity, it 
suffices to show that 

(0A/,0Bt//)=O 

for all / e Ve. Note that 

{QAf,QBUf) = {eAj,QB,Uf) + {eAj,QB,Uf). 

Note also that 



[Uf]{x) = j\f{i)e-'-'^^di 
Jo 



via the variable substitution ^ ^ ^ — 1. 
Therefore, we have that 



= (/(0,-e^%))^^[o,i]- 



Likewise, 



{QAj.QB,Uf) = Y, /'/(e)e-'"*('*+'/')«c?e- /'e2^-/(e)e-2-<cie 
= (e-^^V"(e),/(e))L^[o,i]. 



Therefore, 

Let m denote Lebesgue measure on W^. 

Corollary 1. The samples {Ajz + Pj : j — 1, . . . ,n; z & Z^} are tight on the band E if and 

only if for each a G A \ {0}, either 

1. m{E{^{E + a)) = 0, or 
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2- EjeP, |rfeMj|-ie-2'^*^^-° = 0. 

Proof. Let V and Cj be as in the proof of Theorem^ let Qj = Tp^ipE- The samples {Ajz + jSj : 
j = 1, . . . ,n; z E Z'^} are tight on the band E if and only if equationElin LemmaElis satisfied. 
Note that the condition for a = is automatically satisfied, with K = Yl'j=i I "^^t 4 |~^- For 
a 7^ 0, we must have 

E M)9A^ + «) = E I det A,rie2-^^-«Xi.(0e-'-^^^-«+")Xi.(e + «) 
= 0, 

from which the statement follows. □ 



An immediate corollary of lemma 2 regarding orthogonality is as follows. 

Corollary 2. The samples {AjZ + (3j : j = 1, . . . ,n; z E Z'^} and {Bjz+'jj : j = 1, . . . ,n; z G 
Z'^} are orthogonal on the bands E and F if and only if for all f G Veuf H V 

n „ 

A (better) complete characterization of determining orthogonality of shifted lattices on 
bands E and F is out of reach at this time. Note that in lemma ^ Cp appears in both 
sequences. We wish to consider orthogonal samples with different matrices; the techniques 
used to prove lemma ^ do not readily extend to this case. We conclude, however, with the 
following two special cases: 1) Bj = Aj, and 2) Aj = A and Bj = B. 

Lemma 5. Suppose G and H and ^h,g o,re as in lemma\^ ^h,g = z/ and only if for all 
a e A, 

(6) E I detCp|-^|;;(0^(e + a) = a.e. C 

Proof. If equation IHl holds for all a G A \ {0}, then Qh,g is a Fourier multiplier; if the 
equation also holds for a = 0, then the symbol of flH,G is 0. Conversely, if Qh,g = 0, then 
^H,G commutes with all translation operators Ty, hence equation IHl holds for a G A \ {0}; 
moreover, since flH,G = 0, then the equation must also hold for a = 0. □ 

Let A be as above; for q G M'^, define Iq^a '■= {j £ Va '■ Ij — Pj = q}- 

Theorem 3. The samples A := {AjZ + jSj : i = 1, . . . ,n; z & Z'^} and B := {AjZ + jj : 
j = 1, . . . , ra; z G Z'^} are orthogonal on the bands E and F, respectively, if and only if for all 
a G A either 

1. m{En (F + a)) = 0, or 

2- EjG/, . I det Aj|-ie-2^*^^-" = 0, for all q G R'^. 
In particular, m{E (1 F) = 0. 
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Proof. As above, let V, Cj, and gj be as in Corollary El let hj = T^^ipp- The samples A and 
B are orthogonal on E and F respectively if and only if G^G^ = 0. Let G := {T^ kOj '■ j = 
1, . . . , n; keZ'^} and H := {TA^khj : j = 1, . . . , n; ke Z'^}. Note that fi// g = if and only 
if Q^Ga = 0. 

Applying Lemma |H1 yields that for all a G A, we must have 

J2 I det A,|-'e2-^^^-«Xi^(0e-'^^^^-(«+")xF(e + a) = 0. 

Clearly, if m{E (1 {F + a)) ^ 0, then the trigonometric polynomial 

J2 |detAfc|-^e-2"^^-"e-2"'^« = 0. 

1=1]- Pj k&Iq,a 

The converse follows directly. □ 
For q e M^, define Jq := {j e {1, . . . ,n} : B-^-fj - A-^/Jj = q}. 

Theorem 4. The samples A := {Az + Pj : j = 1, . . . ,n; z E Z'^} and B := {Bz + jj : j = 
l,...,n;z G Z'^} on the bands E and F, respectively, are orthogonal if and only if for all 
m G TJ^, either 

1. m{A*E n {B*F + m)) = 0, or 

2- EjGJ, e-2^^T^-^'"' = 0, for all q G W^. 

In particular, m{A*E n B*F) = 0. 

Proof. The proof follows the proof of Theorem ?? with the following modification. If C 
is a. d X d invertible matrix, then define the unitary dilation operator Dc on L'^{M.'^) by 
Dcg{x) = a/I det C\g{Cx). We have 

n 

Y TAkTp^i)E) {TBkT^^ifjF, g) 

n 

= ^ ^ (-Da/, DATAkTf^^^pE) {DBTBkTy.i^F, Dsg) 
i=i fcezd 

n 

= Y,Y'^DAf.TkDATp^ilJE){TkDBT^,ipF,DBg) 
= {nH,GDAf.DBg) 

where G := {TkDAT^^ipE : j = 1, . . . , n; k e Z'^ and := {TuDbT^j'^Pf : j = I, ■ ■ ■ ,n; k e 
Z'^. It follows, since Da and Db are unitary, that the samples A and i? are orthogonal 
on the bands E and F if and only if ^u,g = 0. Thus, we apply Lemma with A = Z'^, 
gj = DaTp^iPei and hj = DbT^^iPf- The computations are as in Theorem ??. □ 
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